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Abstract 

This is the second in a series of works devoted to small non-selfadjoint 
perturbations of selfadjoint semiclassical pseudodifferential operators in di- 
mension 2. As in our previous work, we consider the case when the classical 
flow of the unperturbed part is periodic. Under the assumption that the flow 
average of the leading perturbation vanishes identically, we show how to ob- 
tain a complete asymptotic description of the individual eigenvalues in certain 
domains in the complex plane, provided that the strength of the perturbation 
e is S> h 1 / 2 , or sometimes only 3> h, and enjoys the upper bound e = 0{h 5 ), 
for some 5 > 0. 
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1 Introduction 

This paper is the second one in a series dealing with the semiclassical spectral asymp- 
totics of small non-selfadjoint perturbations of selfadjoint operators in dimension 2. 
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In our previous work [H^ we studied the case when the classical bicharacteristic flow 
of the unperturbed operator is periodic in some energy shell, and in this paper we 
shall take a next step in this study, still assuming the periodicity of the underlying 
classical flow. 

In this entire project we have to a large extent been motivated by the work of 
A. Melin and the second author where it has been observed that in dimension 
2, there exist wide and stable classes of non- self adjoint operators, for which one can 
obtain a detailed information about the individual eigenvalues in the semiclassical 
limit h — > 0, in some fixed domain in C. While the corresponding conclusion is well- 
known and classical in the selfadjoint case in dimension 1, important difficulties 
appear when passing to higher dimensions, unless one makes an additional strong 
assumption of complete integrability, also on the quantum level. The latter case has 
a long tradition in semiclassical analysis, and we refer to [3] for a recent study of it 
in dimension 2, in the presence of singularities — see also the references given in that 
paper. 

The aforementioned difficulties in the selfadjoint case in higher dimensions are 
intimately related to the phenomenon of exceptional sets, corresponding to unstable 
tori in the KAM theorem, and a basic discovery of jTH] which allowed for complete 
spectral results, was a version of the KAM theorem without small divisors in the 
complex domain. Also in [T^j, an important role was played by certain flow invariant 
Lagrangian tori in the complexified phase space, and we have seen how to obtain 
complete spectral results in some regions of the complex plane, when the strength 
of the non-self adjoint perturbation e = 0(h s ), 5 > 0, is bounded from below by an 
appropriate power of the semiclassical parameter. Corresponding lower bounds will 
play an important role also in this work. 

Our general assumptions will be the same as in jT^j, and as in that paper, we 
shall let M stand for either R 2 or a real- analytic compact manifold of dimension 2. 
We shall also write M to denote a complexification of M, so that M = C 2 when 
M = R 2 , and in the manifold case, M is a Grauert tube of M. 



be the Weyl quantization on R 2 of a symbol P(x, £, e;h) depending smoothly on 
e G neigh (0, R) and taking values in the space of holomorphic functions of (x, £) in 
a tubular neighborhood of R 4 in C 4 , with 



there. Here m is assumed to be an order function on R 4 , in the sense that m > 



When M = R 2 , let 






(1.2) 
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and 

m(X) < C (X - Y) N °m(Y), X, Y e R 4 , (1.3) 
for some C , N > 0. We also assume that 

m > 1. (1.4) 

Assume furthermore that 

oo 
3=0 

in the space of holomorphic functions. We make the ellipticity assumption 

|po, e (z,OI > ^m(Re(x,0), \(x,£)\ > C, (1-6) 

for some C > 0. 

When M is a compact manifold, we let 

Pe= ^Mh){hD x Y (1.7) 

|o|<m 

be a differential operator on M, such that for every choice of local coordinates, 
centered at some point of M, a ate (x; h) is a smooth function of e with values in the 
space of bounded holomorphic functions in a complex neighborhood of x = 0. We 
assume furthermore that 

oo 

a a ,e( x '' h ) ~ ^2 a a,e,j(x)h 3 , h -> 0, (1.8) 
i=o 

in the space of such functions. The semiclassical principal symbol in this case is 
given by 

and we make an ellipticity assumption 

Wx,OI>^<0™ (^,0eT*M, |C|>C, 

for some large C > 0. Here we assume that M has been equipped with some analytic 
Riemannian metric, so that |£| and (£) = (1 + |£| 2 ) 1,/2 are defined. 
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Sometimes, we write p e for p 0t and simply p for p 0>0 . Assume 

P e=0 is formally selfadjoint. (1-10) 

In the case when M is compact, we let the underlying Hilbert space be L 2 (M,/j(dx)), 
where fi(dx) is the Riemannian volume element. 

For h > small enough and when equipped with the domain H(m), the naturally 
defined Sobolev space associated with the weight m (so that in the compact case, 
H(m) is the standard Sobolev space H m with an h dependent norm), P e becomes a 
closed densely defined operator on L 2 (M). Moreover, the assumptions above imply 
that the spectrum of P e in a fixed neighborhood of G C is discrete, when h > 
and e > are sufficiently small. Clearly, if z G neigh(0, C) is an eigenvalue of P £ , 
then Imz = 0(e). As in we shall be interested in the asymptotics of individual 
eigenvalues of P t inside a band \lm.z\ < 0(e). 

Assume for simplicity that (with p = p t= o) 

p _1 (0) n T*M is connected. (1.11) 

Let H p = p^ ■ J| — p' x ■ J| be the Hamilton field of p. In this work we shall always 
assume that for E G neigh (0, R): 

The ifp-flow is periodic on p^(E) n T*M with (1.12) 
period T(E) > depending analytically on E. 

When / is an analytic function defined near p _1 (0) D T*M, we introduce the 
trajectory average, 

T(E) 

f o exp (tH p ) dt on p~\E)nT*M, (1.13) 

which Poisson commutes with p, H p (f) = 0. Let us write out the first few terms in 
a Taylor expansion of p e , 

p f = p + ieq + e 2 r + C(e 3 m), (1-14) 

in the case M = R 2 , and p e = p + ieq + e 2 r + 0(e 3 (£) m ) in the compact case. 
In [Tj)], restricting attention to the range 

/i<e = £>(//), 

for an arbitrary but fixed 8 > 0, we have obtained complete asymptotic expansions 
for all eigenvalues of P e in rectangles of the form 

[-1/0(1), 1/0(1)] + ie[F - 1/0(1), F + 1/0(1)], 
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under the assumption that F is either a regular value of Re (q) , viewed as a function 
on p -1 (0) fl T*M, or a critical value corresponding to a non-degenerate maximum 
or minimum of this function. When the subprincipal symbol of P t= o vanishes, we 
were able to treat even smaller values of e, h 2 <^ e < 0{h s ). 

Now there exist many natural situations when the averaged perturbation (q) 
vanishes identically, and in section 7 of jTH] we have already observed that one such 
case occurs when studying barrier top resonances for the semiclassical Schrodinger 
operator, in the case of a 1 : 1 resonance for the fundamental frequencies at the top 
of the barrier. The purpose of the present work is to carry out a spectral analysis of 
P e under the assumption that (q) = 0. This case has a long tradition in the study 
of self adjoint operators with a periodic classical flow (jSj, jH], [HI], Ell), and at this 
point we would also like to mention an interesting paper by L. Friedlander jHj, which 
studies the spectral localization for the Laplace operator on the n-sphere, perturbed 
by a complex-valued bounded odd potential — see [S] for a corresponding study in 
the real case. A common feature of all the works mentioned above (although not 
of [Hj), which was also exploited in ^Hj, is an averaging procedure, which allows 
one to reduce the dimension by one unit. (For that it is important that we have 
e = 0(h s ), 5 > 0.) Here we also follow this idea, although, as we shall see, some 
modifications in the averaging procedure, as compared with will be required. 
We would also like to emphasize that the implementation of the averaging method 
in the present non-selfadjoint setting involves conjugating our operator by means of 
Fourier integral operators with complex phase, suitably realized on the FBI trans- 
form side, and as in we are naturally led to work in modified exponentially 
weighted spaces, in the spirit of jTH], [13], [TH] . Additional modifications of the 
spaces will be required, and for the final spectral results, we shall sometimes also 
have to introduce a supplementary global hypothesis of dynamical nature, which will 
ultimately allow the spectral problem to become microlocalized to a small neighbor- 
hood of a suitable Lagrangian torus. Let us also remark that we expect that global 
hypotheses of a similar kind will appear naturally in further works in our series, in 
more general situations, when the energy surface of the unperturbed symbol admits 
certain invariant Lagrangian tori with good arithmetic (diophantine) properties. At 
many essential points however, the proofs will be the same as in ^5], and for that 
reason, the presentation of those parts of the proofs will be somewhat less detailed. 

The plan of the paper is as follows. 
In section 2, we reexamine the averaging procedure of [T^j, adapting it to the present 
case, and describe a localization of the spectrum of P e . 

In section 3 we construct a quantum Birkhoff normal form near a suitable torus, 
and, provided that a certain averaged correction has a nondegenerate imaginary 
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part, obtain the first spectral result. This is done completely along the lines of the 
analysis of [15J. 

In section 4, we study the complementary case when the relevant averaged correction 
is real on the real domain, and we then show how to obtain complete spectral results, 
by introducing an additional global dynamical assumption. 

In section 5 we apply the results of section 3 to barrier top resonances, thereby 
complementing the discussion of section 7 of ^3] • 

In section 6, we give a first application of the results of section 4 to the spectrum of 
the damped wave equation on the 2-sphere ( [T], [T7j, [20], [II])- 

Acknowledgment. The first author gratefully acknowledges the support of the 
MSRI postdoctoral fellowship as well as a partial support of the National Science 
Foundation under the grant DMS-0304970. 



2 Averaging reduction and spectral localization 

In this section we shall work near p _1 (0) H T*M and the following arguments will 
be valid in any dimension n > 2. 

Let us recall that in a neighborhood of p~ 1 (0) fl T*M, the principal symbol of 
P e has the form 

p t = p + ieq + e 2 r + 0(e 3 ). 

Assuming that (q) = 0, we shall first reexamine the reduction by averaging, de- 
scribed in section 3 of [12], and adapt it to the present situation. Let Go be an 
analytic function defined near p" 1 ^) fl T*M such that 

H p G = q. (2.1) 

As remarked in [T^j, we may take 

1 fT(E) 

G = —— t(qoexp(tH p ))dt on p~\E)nT*M. (2.2) 

1 {&) Jo 

Replace now T*M by the new IR-manifold A e c = exp (ieHQ )(T*M), which is 
defined in a complex neig hborhood of n T * M - Writing 

(x,£) = exp(ieH Go )(y,7]), 

and using p = (y,rj) as real symplectic coordinates on A eGo , we get 

(Pe)\A tGo = (p + ieq + e 2 r + O(6 3 ))(exp(zeH G0 )(p)) (2.3) 
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k=0 



P + e 2 (r-H Go q- l -H 2 Go p^+0{e 3 ). 



We shall assume in what follows that e = 0(h s ) for some fixed S > 0, and iterating 
the procedure above, we shall modify the choice of the weight function Go and of 
the corresponding IR-deformation, in order to improve ()2.3|) . Set 

G = G + ieG 1 + O(e 2 ), (2.4) 

where G\ is to be chosen, and compute the restriction of p e to the corresponding 
IR-manifold A eG . We get 

p e (exp (ieH G )) = p + ieq + e 2 r + ieH G (p + ieq) - —H G p + 0(e 3 ) 

e 2 

= p + ieq + e 2 r + ie {H Go + ieH Gl ) (p + ieq) - jH GqP + C(e 3 ) 

e 2 

= p + ieq + %e (H Go p + ieH Go q + ieH Gl p) + e 2 r - jH GqP + C(e 3 ) 
= p + e'fr + H^-H^q-^H^+Oie 3 ). 

Using that H p Gq = q, we see that the expression in front of e 2 takes the form 

s ■= r + H p Gi - ^H Go q. 

We choose G\ as an analytic solution in a neighborhood of p _1 (0) PI T*M of the 
equation 

H p G 1 = ^(H ao q-(H aa q))-(r-{r)), 
and then get a reduction of the original symbol to 

p + e 2 (s) + 0(e 3 ), (2.5) 

where 

(*> = (r) ~ \{H Go q) = (r) - {G , q} o exp (tH p ) dt, 
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so that p and (s) are in involution. Recalling the expression for Go from ()2.2|) and 
writing T(p) rather than T(E) for the period on the energy surface p~ 1 (E), we get 



I rT(j>) 

{G ,q} = —— s{qoexp(sH p ),q}ds 
T{p) Jo 

( i rT(p) \ 

+ {T{p), q]\q- Tj^j ^ sq o exp (sH p ) ds I . 

It follows that 

(s) = (r) (2.6) 



2T(p) 
1 

2J\p) 



j--^ J J u{q o exp ((« + v)H p ), q o exp (vH p )} dudv 



Tip) 

q o exp (vH p ){T(p), q o exp (vH p )} dv 



+ 2t\p) j I u ^ T ^ ,q ° exp ( vH p)}Q ° exp (( u + t ') jff p) rfMC ^> 

where the integration in the double integrals in the right hand side of f!2.6|) is per- 
formed over the rectangle [0, T(p)] 2 . Notice also that any other choice of Go satisfy- 
ing (|2.ip gives the same expression for (s). Indeed, if h is a function invariant under 
the Hp-How, then 



I r-i (P) 

({h,q}) = Y(pj J {h,q}oexp(tH p )dt 

I fT(p) 

= — - / {h,qo exp (tH p )} dt = {h, (q)} = 0. 
Tip) Jo 



For future reference, we remark that since (q) = 0, it is true that (Go) = 0. 
Indeed, 

1 rT(E) / rT(E) \ 

(G ) = T( ^ E y J u\j qoexpHu + v)Hp)dvJdu = 0. 

After replacing p e by p e oexp (ieHa), and correspondingly P e by U~ x PJJ ey where 
U e = e eG ( x < hD *)/ h i s a Fourier integral operator quantizing exp [ieHa), which is de- 
fined microlocally near p _1 (0) H T*M, we obtain a reduction of our operator P e to 
an operator with the leading symbol 

Pe =p + e 2 (s) + C(e 3 ). 
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Moreover, if the subprincipal symbol of P e is 0(e), then this is also true after the 
conjugation, in view of the improved Egorov property of U e — see section 2 of |15j . 
As in section 3 of ^5], let g : neigh(0, R) — > R be the analytic function defined 

by 

9 ' {E) = ^r> ^ (0) = ' (2 - 7) 

so that g op has a 27r-periodic Hamilton flow. Set / = g^ 1 . 

Proposition 2.1 Assume that the subprincipal symbol of P t= o vanishes and that 
e <C h 1 ! 2 . Then the spectrum of P e near is contained in the union of the rectangles 
of the form 

4(e) = / fh (k - - A) + [-0(e 2 + h% 0(e 2 + h 2 )] (2.8) 
+ t[-0{e 2 + eh),0{e 2 + eh)}, 

for k G Z. Here a G Z and S G R are the Maslov index and classical action, 
respectively, computed along a closed H p -trajectory C p _1 (0) D T*M. 



Remark. The vanishing of the subprincipal symbol of P e =o in this result could be 
weakened to assuming that the trajectory average of it is constant on each energy 
surface p~~ l {E) n T*M, E G neigh(0,R). Notice also that since e < h l/2 and 
/'(0) > 0, it is true that 4(e) H 4'( e ) = ^> when k ^ k' and /i is small enough. 

Proof: The averaging reduction described above together with the spectral theorem 
shows that it suffices to prove the proposition when e = 0, in which case it is well- 
known and was established in the semiclassical case in ^2], following the earlier 
works [2S1 and |2j. Here we shall take the opportunity to sketch a proof of this result, 
which does not rely upon the averaging procedure on the operator level as in , but 
rather works directly with microlocal representatives near the closed orbits. Similar 
arguments will be used in what follows. In doing so, we shall have to assume that 
T(0) is the minimal period of every closed ifp-trajectory in p _1 (0) (1T*M. It is quite 
likely however, that the following argument can be modified to cover the general case 
when there exist subperiodic orbits, and we plan to return to that in the future. 

Let 7 C p _1 (0) fl T*M be a closed ifp-trajectory of minimal period T(0). From 
section 3 of we recall that there exists a unitary Fourier integral operator, 
microlocally defined from a neighborhood of 7 in T*M to a neighborhood of r = 
x = £ = in T*(Sl x R™ _1 ), such that after the conjugation by this operator, the 
operator P = P e =o takes the form P = f(hD t ) + 0(h 2 ), where the remainder is 



9 



an /i-pseudodifferential operator whose symbol is 0(h 2 ). Here we also use that the 
Fourier integral operator has the improved Egorov property. The operator P acts on 
the space L| of functions defined microlocally near r = x = £ = 0in T*(S\ x R™ -1 ), 
and satisfying the Floquet-Bloch condition 



u(t - 2n, x) = e l( h + 2 4 a) M. 



If now z G neigh(0, R) avoids the union of the intervals 7^(0), then we see that the 
operator 

P-z = f(hD t ) + 0{h 2 ) - z (2.9) 
is invertible, microlocally near r = x = £ = 0, with the norm of the inverse being 

o(h- 2 ). 

Take now finitely many closed trajectories 71, . . . 7^ C fl T*M and small 

open Hp-invariant neighborhoods Qj of jj, 1 < j < N, such that the Q^s form 
a cover of p _1 (0) n T*M. Take also functions Xj e C~(^j)> 1 < j < N , such 
that Y2f=i Xj = 1 near p _1 (0) fl T*M, and the Xj are m involution with p. When 
z G neigh(0, R) avoids the intervals ifc(0), we consider the equation 

(P — z)u = v , u G H{m). 

In what follows we write L^- to denote the Fourier integral operator intertwining P 
and the operator P, microlocally in Qj. Then 

(P - z)UjXjU = Uj (xjv + [P, Xj]u) , 

modulo an 0(h°°)-error. Here [P, Xj] = 0(h 3 ) in the operator sense, since the 
subprincipal symbols of P and Xj both vanish. It follows that, with || • || denoting 
the L 2 -norm, 



Xj u\\<o(^j\\v\\+0(h)\\u 



Summing these estimates over all j and combining them with the usual elliptic 
bound away from p _1 (0) fl T*M, we conclude that 



u\\<0 



1 

V 2 



when h is small enough. The injectivity and hence the invertibility of P — z, for 
z G neigh(0, R) avoiding the rectangles ifc(0), k G Z, follows. □ 
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Remark. If (s) is real on the real domain, then we obtain an improvement on the 
vertical size of the rectangles Jfc(e) in (|2.8|) . and we see that if z G Spec(P e ) fl 
neigh(0, C), then 

kz= 0{e 3 + eh), 

for e = 0(h 5 ), 5 > 0. 

3 Normal forms and eigenvalues in the torus case 
I 

The arguments in this section will to a large extent be parallel to the corresponding 
discussion in sections 4 and 6 of and therefore the following presentation will 
be less detailed. 

We have seen in section 2 that we may reduce ourselves to the case of an operator 
with the leading symbol 

p e =p + e 2 (s) + C(e 3 ). 
Throughout this section, we shall assume that 

Re (s) is an analytic function of Im (s) andp, (3-1) 
in the region of T*M where \p\ < 1/(9(1). 

Introduce the if p -invariant set 

A ,o : p = 0, Im (s) = 0, 

and assume that the minimal period for all the closed ifp-trajectories in Ao,o is T(0), 
and that dp and d(lms) are linearly independent at each point of Ao,o- Then Ao,o is 
a Lagrangian manifold which is also a union of tori. Assume for simplicity that A 0; o 
is connected, so that it is equal to a single Lagrangian torus. In a neighborhood of 
A , p and Im (s) form a completely integrable system, and precisely as in section 

4 of |15j . we perform a real analytic canonical transformation 

k : neigh (f = 0, T*T 2 ) -> neigh (A , , T*M) , (3.2) 

such that p o k = p(£i)> Im (s) o k — Im (s)(£)- (I n section 4 of jT^j we recalled the 
construction of k — see also 5J.) It is therefore clear that when working micro locally 
near A .o and implementing k by means of a microlocally unitary Fourier integral 
operator, which also has the improved Egorov property, we reduce ourselves to a 
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new operator, still denoted by P e , which is microlocally defined near the zero section 
in T*T 2 , and which has the principal symbol 

Also, the transformed operator will act on the space L^(T 2 ) of microlocally defined 
Floquet periodic functions on T 2 , satisfying 

u (x -v) = e i0 - u u(x), v e (2ttZ) 2 , 6=—- + ^-. 

Here S = (Si, S2) with Sj being the action associated to the fundamental cycle 7^ 
in A Q)0 , j = 1, 2, with 71 being given by a closed if p -trajectory, and a = (ati, a 2 ) is 
the corresponding Maslov index. The assumption about the linear independence of 
the differentials implies 

%p(0)^0, %lm(s)(0)^0. 
Using the first of these assumptions and repeating the arguments of section 

TOO 

ju=0 ■ 



4 of [12], we find an /i-pseudodifferential operator A = J2T=oh Ua v(x,£,,t), with 



a = 0(e ), such that the full symbol of the operator 

k 



el A P&-i A = e^P € = ^±(±adA] P e 

1 — n ' \ / 



k=0 

is independent of x\. We remark that the functions a u (x, £, e) are constructed as 
formal power series in e, with coefficients holomorphic in a fixed complex neighbor- 
hood of £ = 0. These formal power series are then realized as C^-symbols since 
we work under the assumption that e = 0(h s ) for some fixed 5 > 0. We are thus 
reduced to the operator 

00 

P e = Y,h%(x 2 ,Z,e), (3.3) 

3=0 

with 

Po=p(^i) + ^(s)(0 + 0(e 3 ), 

also independent of X\. We would like to perform further conjugations of P e by means 
of Fourier integral operators, in order to make its leading symbol independent of x 2 
as well. We shall do so first in the case when the subprincipal symbol of P e= o is not 
necessarily zero. On the symbol level, we write 

Pe = P (^i)+e 2 ((s)(^ + 0(e) + ^pi(x 2 ^,e) + h^p 2 (x 2 ^,e) + ...) (3.4) 

h h 
= p(fi) + £ 2 (r (x2, f , e, -j) + hri(x 2 , £,£,—) + ...), 
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with 

r (* 2 ,£,e, 4) = (*>(£) + 0(e) + %i = <*>(0 + + (4 



( - 1 ' < - ' 1 (- 



Notice that r j = 0{h/e 2 ) for j > 1. 

Following [15] . we shall treat h/e 2 as an independent small parameter. In the 
case when h/e 2 < h Sl , 5\ > 0, an inspection of the arguments of section 4 of [T3] 
shows that there exist operators 



and 



such that 
has a symbol independent of x 



B = b (x 2 , hD x , e, h/e 2 ), b = 0(e + h/e 2 ), 

oo 

Bi = Y, M^2, hD x , e, h/e 2 )h», b u = 0{e + h/e 2 ), 

p e :=e £ad Ble £adB p e ( 3 5 ) 



h, . _ h 



Pe = + e 2 M£, e, ^) + hn{i, e,^) + ...), (3.6) 

with r = (s)(0 + 0(e + h/e 2 ), and r„ = 0(e + h/e 2 ) for u > 1. 

Remaining in the general case, we shall now discuss the construction of the 
conjugating operators assuming merely that h/e 2 < 5 , for some 5 > small 
enough but fixed. Again, this is done completely along the lines of the analysis 
of ^2]- We therefore find that there exists a holomorphic canonical transformation 
K from a complex neighborhood of the zero section in T*T 2 to a similar set, with a 
generating function of the form 

1> \x,r],e, 4) =x-r] + ip per (x 2 ,r],e,^) , ?p pei = O (e + h/e 2 ) , (3.7) 



such that 



(r o K )(y, 7/, e, 9 = (r„(-, £, e, ^)) = (r (-, r/, e, ^)> + G(e 2 + (-) 2 ) (3.i 
e z e A e A e z 

is independent of y. Here the average of r$ is with respect to the x 2 variable. 
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We can quantize Fourier integral operator U and after conjugation by 

this operator, we may assume that we have a new operator P e as in (J3.4)) . with 
r o — ( s )(0 + + V g2 ) independent of x and with rj = 0(e + h/e 2 ), j > 1. 

As before, we can then make a further conjugation e^ ads i in order to remove the 
x-dependence completely, and the conclusion is that if we make no assumption on 
the subprincipal symbol and restrict the attention to h/e 2 < 5q, for <5 > small 
enough, then we can find a Fourier integral operator, 

U-'uix; h) = — j j e^^-y^a(x, rj; h)u(y)dyd V , (3.9) 

with tp{x,r]) = x ■ r] + ifj per (x 2 ,f],e, h/e 2 ), ip per = 0(e + h/e 2 ), and 

°° h h 

B 1 = Y,K(x 2 ,hD x ,e,-)h», b u = 0(e + -), 

such that 

p e ■= e^ ads iAd[/P e 

has a symbol independent of x as in (j3.fij) . with the same estimates as there. 

We shall now turn the attention to the case when the subprincipal symbol of P e= o 
vanishes. The arguments of ^Hl together with the preceding discussion make it clear 
that in this case we may reach the range h/e < 5 , for 5 > small enough. Indeed, 
we write as in ()3.4|) . using that pi = eqi, 

Pe = p{i l ) + e 2 {{s){i) + 0{e) + \ l {x 2 ,ii,e) + %p 2 + h%p, + ...) 

e e z e z 

= p(£i) + £ 2 {r {x 2 , ^, e, -) + hr x {x 2 , l\, e, -) + h 2 r 2 + ...), 



with 



r Q {x 2 ,i,e,-) = (s)(0 + O(e) + -q 1 + f ^p 2 , 
e e e z 

ri{x 2 ,i,e,-) = —p 3 , 



e" e 2 



r 2 {x 2 ,i,e,-) = —Pa,... 



e e- 



It follows that under the stated smallness assumption on e, we can eliminate the x 2 - 
dependence by means of conjugations by Fourier integral operators . To be precise, 
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as in we find that there exist an elliptic Fourier integral operator U 1 and 
Bx(x2,hD x ,e,h/e;h) such that 

eX^AduPt = P(hD X} e, h/e; h) 

has a symbol P(£, e, h/e; h) of the form 

+ e 2 (r + tin + ...), 

where r = (s) + 0(e) + 0(h/e) , and Tj = 0(e + h/e), j > 1. This quantum microfo- 
cal Birkhoff normal form near the torus A 0) o leads to the formal quasi-eigenvalues 
associated with A ,o and given by P(h(k — a/4) — S/2ir, e,h/e ; h), k £ Z 2 . 

So far our discussion has been completely analogous to the corresponding one in 
section 4 of ^H], and it is also clear that the analysis of the globally well-posed 
Grushin problem from section 6 of that paper can be applied without any change. 
To fix the ideas, in what follows, let us discuss the spectral results assuming that the 
subprincipal symbol of P €=0 vanishes and h/e < 8 <C 1, e = 0(h 5 ), 5 > 0. We then 
inspect all the steps of the argument of section 6 of [THj to find a new (/i-dependent) 
globally defined Hilbert space H(A e ), associated with a suitable IR-manifold A e , 
(e + /i/e)-close to T*M, such that the action of P t on H(A e ) can be, microlocally 
near A 0) o, identified with the action of P(hD x , e, h/e; h) on L^(T 2 ), by means of an 
elliptic semiclassical Fourier integral operator 

U = 0(1) : H(A e ) - L 2 (T 2 ). 

We refer to Proposition 6.1 of ^3] for the precise description of the basic properties 
of U. As in , the construction of the deformation A e is done in such a way that 
for p G A e away from a small neighborhood of Ao,o, 

\ReP e (p,h)\>-^- or llmP^)^^-. (3.10) 

Let now z vary in the rectangle 







|Rez| < — , 









for a sufficiently large constant C > 0. If z in (jH.llj) is away from any e 2 h/C- 
neighborhood of the values P(h(k — a/4) — S/2tt, e, h/e; h), k e Z 2 , it follows from 
the argument explained in section 6 of that the operator 

P e - z : H(A e ) -> H(A £ ) 
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is invertible, with the norm of the inverse 0((t 2 h)^ 1 ). The setup of the Grushin 
problem is also identical to the one in jT^j, and repeating the arguments of that paper 
we find that up to 0(h°°), the eigenvalues of P t in (J3.ll)) are given by the quasi- 
eigenvalues associated to A ,o, P(h(k — f ) — j-, e, -; h), k e Z 2 . The corresponding 
result is also true in the range h 1 / 2 < e = 0(h s ), 5 > 0, if the subprincipal symbol 
of P e =o is not necessarily zero. 

We summarize the results of this section in the following theorem. 

Theorem 3.1 Let P e be an operator as in the introduction, so that (j!.12j) . 1)1.14)1 
hold. Assume that the trajectory average (q) vanishes. Define s as in section 2, so 
that (s) is given by (12.6)) . and assume 1)3.1)) . Let Fq be a regular value oflm(s), 
considered as a function on p _1 (0) D T*M. Assume that the Lagrangian manifold 

A ,F : P = 0, Im (s) = F 

is connected, so that it is diffeomorphic to a single torus, and that T(0) is the 
minimal period for the H p -flow in Ao,f - We fix a basis for the first homology group 
o/Ao,f given by the cycles jj, j = 1,2, with 71 corresponding to a closed H p -orbit. 
Let us write S = (Si, S2) and a = (at, 0^2) for the actions and Maslov indices of the 
cycles, respectively. Assume that h 1 / 2 C e = 0(h s ) for some 5 > 0, and let C > 
be large enough. A complex number z satisfying 







\Rez\ < —, 


Im z\ 







is in the spectrum of P e precisely when 

z = z (k, e;h) = P (h { k -j)-^^^ h )+ °( h °°)' k e 7? - 

The function P(^,e,^;h) is holomorphic in £ e neigh(0, C 2 ), smooth in e, 4r G 
neigh(0,R) ; and as h — > 0, has an asymptotic expansion in the space of such func- 
tions, 

P(C, e,^;h)~ p(e0 + e 2 (r ^, e, £j + h n ^, e, + . . .) . 

We have 

r =(8K£) + o(e + £), rj = o(e+^y j > 1. 

In the case when the subprincipal symbol of P t=0 vanishes, we have a correspon- 
ding result for e in the range ft < e = 0(h 5 ), 5 > 0, with the only difference that 
"h/e 2 " in the description of the eigenvalues should be replaced by "h/e". 
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4 Spectral asymptotics in the torus case II 

Throughout this section, it will be assumed that the subprincipal symbol of P e=0 
vanishes and that (q) = 0. 

Recall that the principal symbol of the original operator P e has the form 

p e =p + i eq + e 2 r + ie 3 w + 0(e*), (4.1) 

in a complex neig hborhood of p _1 (0) n T*M. In some of the applications that we 
have in mind (see section 6), the leading perturbation q is real on the real domain, 
and in this section we shall assume in addition that the next order term, r, is real 
as well. The second averaged correction (s), introduced in (|2.fi|) . is then purely real, 
and the results of the preceding section do not apply. The purpose of this section is 
to understand what kind of spectral results it is possible to obtain then. As before, 
our starting point will be the averaging method. To be precise, we shall now have 
to modify further the choice of the IR-deformation, given by the weight G in (|2.4|) . 
Let us consider 

G = G + ieG 1 + e 2 G 2 + O(e 3 ), (4.2) 

where G 2 is to be determined, and as before, compute the restriction of p t to the 
corresponding manifold A eG = exp {ieH G ){T* M). We get 

(p + ieq + e 2 r + ie 3 w)(exp {ieH G )) = p + ieq + e 2 r + ie 3 w + ieH G (p + ieq + e 2 r) 

(ie) 2 (ie) 3 
+^-H G (p + ieq) + -—J—H G p + C(e 4 ) = p + ieq + e 2 r + ie 3 w 
2 o 

+ie (H Go p + ieH Go q + e 2 H Go r + ieH Gl p - e 2 H Gl q + e 2 H G2 p) 



V 



-- (H 2 GqP + ieH Go q + ie(H Go H Gl + H Gl H Go )p) + ^f-Hfo + 0(e 4 
+ e 2 (-H Go q - H Gl p + r - \h 2 g ^ 

+k 3 (w + H Go r - H Gl q + H G2 p - -fl^p - ^H 2 Go q - ^(H Go H Gl + H Gl H Go )p^j 

As before, using that H p Gq = q, we see that the expression in front of e 2 takes the 
form | 

r + H p Gi - -H Go q, 
and the coefficient in front of ie 3 becomes 

w - \ H G^q - H P G 2 - 2 H G Q - \h Go h Gi p + H Go r. 
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From section 2 we then recall that we choose G% as an analytic solution in a neigh- 
borhood of p _1 (0) fl T*M of the equation 

H p G 1 = ~(H Go q-(H Go q))-r+(r), 

Notice also that G± is real-valued on the real domain, since we assumed that this is 
the case for q and r. Then we see that the coefficient in front of ie 3 is of the form 

t:=w- ^H Gl q - H P G 2 - —H Go q - ^H Go (H Go q) + hl G j + -H Go (r), 

and it is clear that we can choose G 2 , defined near p -1 (0) fl T*M, such that after 
composing the leading symbol p e , which has the form (|4.1|) . with the complex cano- 
nical transformation exp (ieH G ), with G of the form (|4.2j) . we get a reduction of p e 
to 

p + e 2 (s)+ie 3 (t) + 0(e 4 ), (4.3) 

with (s) and (t) real-valued on the real domain, if we also assume that w has this 
property. We shall now give explicit expressions for these corrections, assuming for 
simplicity that the period for the f/p-flow, T(E) = T is independent of the energy 
E e neigh(0,R). The expression (J2.6|) for (s) from section 2 then simplifies and we 

get 

l.„ . , , 1 rT 



(*) = ( r ) - -ji H G q) = (r) - 2p ^ {Co, q] ° exp (uH p ) du, 
and recalling the expression for Go we get 

( s ) = ( r )__L J J u{qoexp((v + u)H p ),qoexp(vH p )}dvdu, (4.4) 

where the integration is performed over the rectangle [0,T] 2 . We also get 

<*> = W-~(H Gl q)-^(H Go q)-\(H Go (H Go q)) 
+ l(H Go r) + l(H Go (r)), 



2 



where 



[ t ((H Go q — 2r) o exp (tH p )) dt on p~ 1 (E)nT*M. 
Jo 
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We shall now simplify the expression for (t). First of all, notice that 

(H Go {H Go q)) = 1 J\g o exp (tH p ), (H Go q)}dt = {(G ), (H Go q)} = 0, 

since we have already observed that the flow average of Go vanishes. Similarly, 

(H Go (r)) = 0. 

If we introduce 

/ = ~Y^2 J J uv{qoexp(vH p ),{qoexp(uH p ),q}}dudv, (4.5) 
g = -Jj^ J J uv{{qoexp((u + v)H p ),qoexp(uH p )},q}dudv, (4.6) 
1 f T 

k= 2TJ M ({9°exp(wi/p),r} + {roexp(w# p ),g}) du i ( 4 - 7 ) 
(t) = (f + g + w + k). (4.8) 



and 



then 



In what follows, we shall consider, as we may, the operator P e , defined microlo- 
cally near p _1 (0) fl T*M, with the principal symbol of the form 

p + e 2 (s)+te 3 (t} + 0(e 4 ), 

with (s) and (t) real- valued, and with the subprincipal symbol 0(e). When deriving 
spectral asymptotics for P t , we shall distinguish the following cases, depending on 
the size of e: 

1. /i < e < h 1 / 2 . 

2. h 1 / 2 < e = 0{h s ), 5>0. 

3. e ~ h 1 ' 2 . 

1. T/ie case when /i<e<C /i 1 / 2 . 

In this case we know, thanks to Proposition 2.1, that the spectrum of P £ near 
has a cluster structure. 
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Let us fix the energy level (0, 0) and introduce the set 

A ,o :p = 0, (s) = 0, 

which is invariant under the if p -flow. We assume that dp and d(s) are linearly 
independent at every point of A , and then, assuming that A is connected, we see 
that this set is a Lagrangian torus. We shall also assume that T(0) is the minimal 
period of every closed if p -trajectory in A 0) o- In a neighborhood of this torus, p 
and (s) form a completely integrable system, and we then pass to the action-angle 
variables by means of a real analytic canonical transformation 

k : neigh (f = 0, T*T 2 ) -> neigh (A , , T*M) , 

such that p o k = (s) o k — (s)(£)- Notice that when expressed in terms of 

the action coordinate £i, p becomes = where the function / has been 

defined after ([2.7)1 . Implementing k by means of a microlocally unitary Fourier 
integral operator, we get a new operator, still denoted by P e , microlocally defined 
near the zero section in T*T 2 , which has the leading symbol 

P(6) + e 2 (3)(0+^ 3 (t)(x 2 ,0 + O(6 4 ). 

The assumption about the linear independence of the differentials implies that 

<%P(0)^0, %( S )(0)^0, 

and repeating the arguments of section 3, we then find an /z-pseudodifferential op- 
erator A = Yl^Lo h u a u (x, £, e), with a = C(e 4 ), such that the full symbol of the 
operator 

-A : 



OO r . \ K 

'P t e-i A = e^ dA P e = J2(T^A) P t 



e h 

is independent of %\. This is done by solving transport equations in x\, obtained by 
looking for each a u (x,£,e) as a formal power series in e. This reduces the following 
discussion to the operator 



oo 



j=0 

with 



Po = p(fi) + e 2 (^)(0 + ze 3 (t>(* 2 , + 0(e 4 ) (4.9) 

also independent of xi, and p\ = 0(e). It is then clear that under the assumption 
that 

h 

-<h Sl , 5i>0, 
e 
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by repeating the arguments of section 3, we can perform further conjugations of P e 
by means of Fourier integral operators, and make its full symbol independent of x 2 , 
modulo 0(h°°). To be precise, we find that there exist 

B = b (x 2 ,hD x ,e,^;h) b = O (e + - 



and 



such that 



h ( h 

B 1 = J2 hU bu(x 2 ,hD x ,e,-), b„ = OU+- 



es ad Bi e £ ad Bo p 



has the full symbol independent of x, modulo 0(h°°). In the case when we only 
have ft/e < 1, as before, the operator e tB °^ h has to be replaced by a more general 
Fourier integral operator, constructed as in section 3 and as in [T3]. We then have 
the associated quasi-eigenvalues given by 

pfh(k-e),e,^;hY fc6Z 2 , # = _L + «. 

Analyzing this microlocal reduction to a normal form near Ao,o, similarly to what 
we did in section 3 and in section 6 of |15j, we get the following result. 

Proposition 4.1 There exists an IR-manifold A C T*M, which is (e + h/e)-close 
to T*M and agrees with T*M away from p _1 (0) H T*M, and a smooth Lagrangian 
torus A 0) o C A ; such that for p G A in a neighborhood of A 0)0 , we have 

p = exp (ieHa) k k{x, £), (x, e T*T 2 , 

and A 0j o corresponds precisely to the zero section £ = 0. Here k is the canonical 
transformation given by the action-angle variables, and k is constructed as before, 
using the generating function (J3.7J) . with h/e 2 replaced there by h/e. When p G A 
is away from a small neighborhood of Ao,o in A and |Re P e {p] h)\ < 1/C for a large 
enough C , it is true that 

\{»){p)\> 1 



0(1) 

Furthermore, there exists an elliptic uniformly bounded Fourier integral operator 
U : H(A) — > Lg(T 2 ) such that, microlocally near Ao,o, UP e = PU . Here 

hD x ,e,^h\ : L 2 (T 2 ) - L 2 (T 2 ) 
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has the full symbol holomorphic in £ ~ in C 2 , and depending smoothly on e, 
h/e G neigh(0, R), 

e, fc) = p(fc) + e 2 (V (£, e, *) + e, + . . . 

r =(s)(Z)+o(e+*Y rj = o(e+*Y j > 1. 

Before setting up a suitable globally well-posed Grushin problem in the Hilbert 
space if (A) in order to identify the spectrum precisely, following [TH] . we shall first 
show that when 2 G neigh(0, C) is such that 



Re,-/(M^-f)-| 



6 2 



(4.10) 



for some k\ G Z, and z avoids the union of the pairwise disjoint open discs Dk 2 (h) 
of radii e 2 h/0(l), that are centered at the quasi-eigenvalues 

\ 4 27T 4 27T e 

for fc 2 G Z, then the operator 

P £ - z : if (A) -> if (A) 

is bijective. In doing so, we shall essentially follow an argument from section 6 
of [Ej, with a slight improvement. Let % £ QT(A) be supported in a small flow 
invariant neighborhood of A 0j o where P e is intertwined with P, and x = 1 near A 0i o- 
We also assume, as we may, that on the operator lever, 

[P e , X ] = 0(h°°) :ff(A)-£T(A). (4.11) 

Indeed, to achieve the property (j4.11|) it suffices to choose a function xo £ C^°(T*T 2 ) 
with suppxo close to £ = and xo = 1 in a neighborhood of this set, such that 
Xo — Xo(0 depends on £ only. Conjugating xo by means of the microlocal inverse V 
of U and modifying it by an (9(/i°°)-factor outside a neighborhood of A 0i o, we obtain 
X with the desired properties. 

Let us introduce a partition of unity on A, 

1 = X + ^1,+ + ^1- + ^2,+ + T/>2,— 
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Here the functions ipi t ± G C£°(A) are supported in flow invariant regions Q±, such 
that ±(s) > 1/(9(1) in Q±, respectively. Moreover, we can arrange so that ipi.± are 
in involution with p, the principal symbol of P e =o on if (A). Finally, ^2,± are such 
that ±ReP e > 1/(9(1) in the support of ip2,±- We shall prove that 

\\(l- x )u\\<o(^\ \\v\\ + 0(h°°)\\u\\. (4.12) 

Here the norms are taken in H(A), (P e — z)u = v, and z G C satisfies f!4.10|) and 
avoids the discs Dk 2 (h), k 2 G Z. When establishing (|4.12|) . we shall first prove this 
bound with in place of 1 — X- 

When A" G N, let ip -< ipi -< . . . -< ip^, ip = tpi t +, be smooth cutoff functions 
supported in Q + , which Poisson commute with p. (The standard notation f ~< g 
means that g — 1 in a neighborhood of supp/.) We have 

(P e - z)^ju = if/jv + [P e , iij}u, < j < N. 

Now the operator P t — z is invertible, microlocally in fi + , with the norm of the 
microlocal inverse being 0(e~ 2 ). It follows that 

\\^u\\<o(^j (\\v\\ + || [^,^11) + 0(h°°)\\u\\. 

Notice that here [P t ,ipj] = 0(h 3 ) + 0(e 2 h) = 0(e 2 h), since the subprincipal symbol 
of P e=0 vanishes and h < e. Now ^(1 — V'i+i) — 0(h°°) in the operator sense, and 
hence we obtain 



[P e , ^j]u 1 1 < 0{t 2 h) 1 1 i/) j+1 u 1 1 + 0{h°°) 1 1 u 



and therefore, 



Combining these estimates for j — 0, . . . N, we get 

1 1 1 1 < O J 1 1 v 1 1 + 0^(1)^1 1 4> N u 1 1 + 0{h°°) \\u\\. 

The same estimate is then obtained for ipi-U, which is concentrated in the region 
where (s) < —1/(9(1), and a fortiori such estimates also hold in regions where 
±ReP e ~ 1. The estimate follows. 
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Relying upon ([4.12)1 . we shall complete the proof of the fact that the spectrum 
of P e in the set (|4.1U|) is contained in the union of the discs Dk 2 (h). Let us write 



(Pe ~ Z)XU = XV+ [Pe, X]U, 

and from (j4.11|) we recall that the norm of the commutator term does not exceed 

0(h°°)\\u\\. 

Applying the Fourier integral operator U of Proposition 4.1, we get, modulo 0(h°°)- 
errors, 

[P-zjUxu = U O + [P e , x)u) ■ 

Now an expansion in Fourier series shows that the operator P — z is invertible, 
microlocally near £ = 0, with a microlocal inverse of the norm (9(e _2 /i _1 ), provided 
that z in the set (|4.1()jl avoids the discs Dk 2 {h). We get 

and combining this estimate together with ()4.12|) we infer that the operator P e — z : 
H(A) — > H(A) is injective, hence bijective, since it is a Fredholm operator of index 
zero by general arguments. Now H(A) agrees with L 2 (M) as a space, and we 
conclude that z is not in the spectrum of P e . 

When z varies in the disc Dk 2 (h) contained in the set (|4.10J) . for some k2 G Z, we 
shall recall briefly the setup of the global Grushin problem for the operator P e — z. 
Introduce the operators R + : H(A) — > C and i?_ : C — > H(A) by 

R + u — (Uxu\ek), R-.U- — uJVe^ fc = (fci,fc 2 )- 

Here 

e k = ^exp (l-(h(k - \ ) - |)arj , 

the scalar product in the definition of R + is taken in Lg(T 2 ), and V is the microlocal 
inverse of U. It is then clear that the arguments of section 6 of apply as they 
stand and show that for every (v,v+) G H(A) x C, the Grushin problem 

(P e — z)u + = v, R + u = v + 

has a unique solution (u,u_) G H(A) x C given by 

u = Ev + E + v + , U- = E-V + E^ + v + . 
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Here the zeros of E h (z) agree with the eigenvalues of P e in D^ 2 (h), and repeating 

the arguments of [15J we find that E-+(z) = z — P(h(k — 9i),h(k 2 — 62), e, j]h) 
modulo an error term which is 0(h°°). We get the following result. 

Theorem 4.2 Let P e be an operator as in the introduction, so that ()1.12|) . (j4.1j) hold. 
Assume that q, r, w are real-valued and that the trajectory average (q) vanishes. 
Define s by (12. 6|) and recall that p e can be reduced by averaging to (|4.3j) . where (t) 
is given by (|4.8J) . (14. 5 J) - (14. 7j) . in the case when T(E) = T is independent of E. 
Assume that the subprincipal symbol of P e= o vanishes. Let Fo be a regular value of 
(s), considered as a function on p (0)f\T*M. Assume that the Lagrangian manifold 

A ,F 'P = 0, (s) = F 

is connected and that T(0) is a minimal period of every closed H p -trajectory in A ,f • 
When 71 and 72 are the fundamental cycles in A 0j f with 71 being given by a closed 
H p -trajectory, we write S = (Si, £2) and a = (0:1,0:2) for the actions and Maslov 
indices of the cycles, respectively. Assume that 

Let C > be large enough. Then for each k\ G Z, the eigenvalues of P e in the set 



Re.-/(^i-f)"§)-^o 



are given by 



pfhih-^-^hih-^-^e^h^+Oih^), hEZ. (4.13) 

Here P(£, e, h/e; h) is holomorphic in £ G neigh(0, C 2 ), smooth in e, - G neigh(0, R), 
and has a complete asymptotic expansion in the space of such functions, as h — > 0, 

P(t,e, -; h] ~ /(£i) + e 2 ( r e, -) + hr x (^e, h 

We have 



ro(O = <«>(O + rj = 0{e+*), j>l. 
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Remark. It is understood that in ()4.13|) we only consider those values of k% G Z for 
which the first argument (£1) of P is small enough. 

Remark. The arguments presented in this case can also be applied to a class of 
self adjoint perturbations of P e =o, with the leading symbol of the form p+eq+O (e 2 m) , 
and with the subprincipal symbol 0(em). (In the compact case, the weight m should 
be replaced by (£) m -) Here q is real on T*M and we assume that the flow average 
(q) does not vanish identically. If e -C h then we know from Proposition 2.1 that the 
spectrum of P e in neigh(0, C) has a cluster structure. Let us restrict e further and 
assume that h 2 <C e. If Fq is a regular value of (q) along p _1 (0), then it follows from 
the arguments of this subsection that we may describe precisely the asymptotics of 
the individual eigenvalues of P e inside the sub-clusters given by 

< C>1, k\ G Z. 

Notice, in particular, that if M is a smooth Zoll surface, and P = —A + XV is a 
Schrodinger operator on M, then the operator h 2 P is of the form above, provided 
that the coupling constant A in front of the potential is sufficiently large, but fixed. 
Also, notice that in this case, the analyticity assumptions are not needed, as we 
may work with Fourier integral operators with real phase and hence stay within the 
framework of the usual L 2 -spaces. To the best of our knowledge, such results in the 
selfadjoint case do not seem to have been formulated explicitly in the literature — see, 
however, [22] for asymptotics of pair correlation functions on Zoll surfaces. 

2. The case when h 1 ' 2 < e = 0(h s ), 5>0. 

Let us recall that by means of an averaging procedure, we have reduced the 
original operator P e to an operator with the principal symbol 

p + e 2 (s)+ie 3 (t) + (D(e 4 ), (4.14) 

where (s) and (t) are real on the real domain. Moreover, the subprincipal symbol 
of the averaged operator is 0(e). 

In the first case, we saw that due to the cluster structure of the spectrum, the 
imaginary part of the principal symbol played no major role in the arguments. The 
situation will be quite different now. When T > 0, we introduce the double average, 

({t)) T =^J (t) o exp (uH {s) ) alu, (4.15) 

defined in a neighborhood of p _1 (0) D T*M. Let G = Gt be a real- valued analytic 
function defined near p _1 (0) H T*M, such that 

H,,)G = (t) - ((t)) T . (4.16) 



H«-/K*.-i)-i)-tfi 
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The equation ()4.16|) is a convolution equation along the if( s }-trajectories, and we 
solve it by putting 

G = J k(u/T)(t) o exp (uH( s ))du, 

where A; is a piecewise linear function, vanishing outside a bounded interval, and 
solving k' = — 5 + l[o,i]- Notice that H p G = 0. After composing the principal 
symbol (j4.14|) with the holomorphic canonical transformation exp(ieHc), and, on 
the operator level, after a conjugation by means of the corresponding microlocally 
unitary Fourier integral operator which also has the improved Egorov property, we 
may assume that our operator is microlocally defined near p _1 (0) n T*M and has 
the principal symbol 

p + e 2 (s)+te 3 ((t)) T + T (e 4 ). (4.17) 

We also know that the subprincipal symbol of the new operator is Ot{c), and we 
can be more precise and remark that it is in fact 0(e) + (9r(e 2 ). 

We keep the basic assumption of case 1 that dp and d(s) are linearly independent 
along the set Ao,o : P = 0, (s) = 0, which is then a flow-invariant Lagrangian torus, 
provided that it is connected. Assume also that T(0) is the minimal period for the 
Hp-How in Ao,o- Consider now the restriction of ((£))t to A ,o- Then the limit 

((*»«, := lim ((t)) T (p), pGAo.o, (4.18) 

T — >oo 

exists, ((t))oo ^ R-j an d the convergence is uniform on A 0j o- We also remark that 
{{t))oo is precisely the mean value of (t) over A 0j o and that 

((t)) T (p) ~ «f»oo = O f^j , (4.19) 

uniformly on Ao,o- Let us also recall the Lagrangian foliation given by the flow- 
invariant tori 

A E:F :p = E } (s)= F, 
for (E, F) E neigh(0, R 2 ). We shall assume that for < a 1, the set 

p- 1 (0)\ |J A , F 

\F\<a 

splits into two disjoint flow-invariant components, which we denote by A a and A_ a . 
We then introduce the following global assumption: 

For any b G (0, a) there exists (4.20) 
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T(b) > and C(b) > such that 

inf («t» T (p) - «*»„) > T > T(b), 

sup {({t)) T (p) - < T > T(b). 



A[- a ,-b] 



|J A ,f]U A - 

<F<a / 



Here we have put 

A[6 )0 ] 

\b<F<a 

and A[_ 0) _6] is defined similarly. 

Given an operator with the principal symbol (|4.17jl . considered microlocally near 
Ao,o, we pass to the standard torus T 2 and eliminate the Xi-variable, exactly as 
before. We are then reduced to the operator P e = P e (T), acting on L^(T 2 ), which 
is defined microlocally near £ = in T*T 2 , and which has the form 

oo 

P e = £V&(a; 2 ,£,e,T), (4.21) 

j=0 

with 

p ( X2 ^^,T)=p^ 1 ) + e 2 (s)(0 + ^((t))T(x 2 ,0 + O T (e 4 ), 

and 

Pi(x 2 ,£,e,T) = egi(ac 2 ,f,e,T) = egi(ac 2 ,f,0) + T (e 2 ). 

Let us write 

Pe = + e 2 (r (x 2 , £, e, p T) + /in(x 2 , £, e, ^ T) + . . .J , (4.22) 

where 



= (s) (0 + ie((t)) T (x 2 , + C?r(e ) + ~gi(x 2 , £, 0) + -p 2 (x 2 , £, e, T), 



and 

n = O t (i) + — p 3 , 



r 2 = 

so that j-j = Ct(1), j > 1. 



r 2 = -rp 4 , 
e 2 



28 



We shall reexamine the construction of the operator B = b (x 2 , hD X) e, ^,T) 
such that 

r oexp (ifa,, ) = 2^ 



fc! 

fc=0 



is independent of x 2 , modulo 0(h°°). As before, we are looking for 6 i n terms of a 
formal power series in h/e — 0(h 1 / 2 ) and e = 0(h 5 ), 5 > 0. Writing 

7 7 7 h , n ( , h. 

% = o ,i, e + 60,0,1- + ( e > 



we find that 6 ,o,i should solve the transport equation 

d&{s)d X2 60,0,1 = 5i(»2,C,0) - (£[(•, £,0)), 

with the average of q± in the right hand side standing for the average with respect 
to the x 2 variable. Therefore, 6 ,o,i = 0(1). As for the term 60,1,0; we see that it 
satisfies 

^<a>9x a 6d,i,o = <(««»T(x2,0-<«*»T», 

where again the last average in the right hand side is with respect to x 2 . Since the 
right hand side here is 0{1/T), uniformly in x 2 and £ G neigh(0, R 2 ), we get 



60,1,0(^2,0 = 



In a similar way, we construct the lower order terms, and obtain 

It follows that there exists 

B = b (x 2 , hD x , e, T), 6 = O (1) + (7) + t ((e, ^) 2 ) , 

and 

5i = e ' -' T )' h » = ° t ^> 



such that 



e 



p _ e 7T ad B le £ad So 



29 



has the full symbol independent of x, 

P = P(£0 + e 2 (r (c, e, 7, + fcn , e, 7, Tj + . . . 

with 

ro = <*>(£) + ie<«*»T>(fl + ° (7) + ^ K + 2 

= Ot(1), J>1- 

Repeating the arguments of section 6 of [T5], we now introduce a globally defined 
IR-manifold A C T*T 2 , with T 2 standing for the standard complexification of T 2 , 
such that A = exp (F 6o )(neigh (f = 0, T^T 2 )) near £ = in T*T 2 , A = T*T 2 further 
away from £ = 0, and such that along A, 

Im6 = 0, Imf 2 = O j + O t (( £ ,M, (4.23) 

Im, = O (l) +O 0) +OT (( f ,^). 

Then the action of P e on the space H(A), associated to A by means of the FBI- 
Bargmann transformation on T 2 , is microlocally near £ = 0, unitarily equivalent to 
the conjugated operator which acts on L 2 (T 2 ) and whose Weyl symbol is indepen- 
dent of x and has the form 

p(&) + e 2 (r (£, e, 7, T) + hn(£, e, T) + . . . 

with the same estimates on rj as above. 

We consider the imaginary part of P e along A, 

ImP £ = e 3 ((t)) T (Rex 2 , Ref) + O (0 + O(efc) + r (e 4 ). 

Choosing now first T sufficiently large but fixed, and then h small enough, and using 
the fact that h 1 / 2 <C e together with the global assumption (|4.20|) . we see that away 
from a small neighborhood of the torus, the absolute value of lmP t — e 3 ((t)) OQ is 
bounded from below by e 3 /(9(l), provided that ReP e < 1/0(1). 

Corresponding to_Jhe manifold A on the torus side, we get a globally defined 
IR-manifold A C T*M, which is an (e + /i/e)-perturbation of T*M, and near A 0; o it 
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is obtained by replacing k(T*T 2 ) there by k(A). Here k is the action-angle canonical 
transformation, and for simplicity here we did not include the complex canonical 
transformations exp (ieHa), coming from the averaging procedures along the flows 
of Hp and . 

We summarize the discussion above in the following proposition. 



Proposition 4.3 Assume that the subprincipal symbol of P e =o vanishes, and that 
the global assumption (|4.2U|) holds. Let e satisfy h 1 / 2 < f = 0(h s ), 6 > 0. Then 
there exists an IR-manifold A C T*M and a smooth Lagrangian torus A 0; o C A 
such that when p E A is away from a small enough neighborhood of Ao,o in A and 
| Re P e {p, h)\ < 1/C for C > large enough, it is true that 



Im P e (p, h) 



«*»c 



> 



0(1) 



(4.24) 



The neighborhood can be taken as small as we wish provided that T in f!4.20|) and 
the implicit constant in (J4.24j) are chosen large enough, and h is sufficiently small. 
The manifold A is an (e + h / e) -perturbation ofT*M and agrees with this set away 
from p _1 (0) R T*M. Microlocally, near Ao,o? the operator P e : H(A) — > H(A) is 
conjugated by means of an elliptic uniformly bounded Fourier integral operator to 
the translation invariant operator P on L 2 e (T 2 ). 



Using this proposition and the methods of section 6 of [T3], we get the final 
spectral result of this case. 



Theorem 4.4 We keep the general assumptions of Theorem 4.2 and consider the 
range h l l 2 C e = 0(h s ), 5 > 0. Assume that the averaged corrections (s) and 
(t) are real on the real domain, and that the differentials of p and (s) are linearly 
independent along the set 

A , :p = 0, (s)= 0, 

which is assumed to be connected. We assume that T(0) is the minimal period of 
every closed H p -trajectory in Ao,o, and we write S E R 2 and a E 7? to denote the 
actions and Maslov indices along the standard fundamental cycles in A 0i0 . Introduce 
((t))oo as the mean value of (t) over A 0i o and make the global assumption (J4.2(Jj) . Let 
C > be large enough. Then the eigenvalues of P e in the rectangle 

1 e 3 
|Rez|<-, \lmz-e 3 ((t))J<- 
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are given, modulo 0(h°°), by the quasi- eigenvalues 

Here P(£, e,-;h) is holomorphic in £ G neigh(0, C 2 ), smooth in e, - G neigh(0, R), 
and has an expansion as h — > 0, 

P(e, e, -; fc) = p(6) + e 2 (r (£, e, -) + hr^, €,-) + ... 
e \ e e 

We have 

r (t,e,*) = (8)($+o(e+!?\, ^ = 0(1), J > 1. 



3. T/ie case when e ~ /i 1 / 2 . 

Assume, as we may, that we are now in a situation when the first two averaging 
procedures have been carried out, so that we are dealing with the operator P t of 
principal symbol 

p , e =p + e 2 (s}+te 3 ((t}} T + O T (e*) } (4.25) 
and subprincipal symbol 

p lj£ = eg 1 + C T (e 2 ), 

where gi does not depend on e and T. Here we may also assume that the subprincipal 
symbol has already been averaged along the flow of p, to the leading order, so that 
Qi = (li)- Since in the relevant parameter range, we have e 3 ~ eh, we shall have to 
consider not only long-time averages ((t))r, but we shall also average (gx) along the 
if( s )-flow. Let therefore G solve 

H {S) G = ( qi ) - (( qi )) T , 

where 




and compose P e with exp (JlHg)- Then, using also the fact that G Poisson commutes 
with p, and implementing exp {-He) by means of a Fourier integral operator with 
the improved Egorov property, we get a new operator with principal symbol ()4.25|) 
and with subprincipal symbol 

e(( qi )) T + T (e 2 ). 
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In what follows, we assume that (0, 0) is a regular value of the mapping T*M 3 
P l— * {p{p)i ( s )(p)) £ R- 2 ) with the connected pre-image, and we shall work microlo- 
cally near the Lagrangian torus 

A ,o : V = 0, (s) = 0. 

As before, we shall also assume that T(0) is the minimal period of the Hp—Row in 
A ,o- 

Introducing the quantities ((t))oo an d ((<?i))oo, defined as the mean values over 
Aq,o of (t) and (qi), respectively, we recall next that as T — > oo, 



((t))T = ((t))oo + ((ll))T = «gi»oo + Ol- 

uniformly along Ao.o- Using the notation of case 2, we introduce now the following 
global assumption, which is a direct analogue of the hypothesis (|4.2()jl : 

For any b £ (0, a) there exists (4.26) 
T(6) > and C(b) > such that 

inf (((lm qi )) T (p) - ((lmq 1 )) QO ) > T > T(b), 

A[6, a] V / 6(D) 

^sup (((Im gi )) r (p) - ((lm qi ))Jj < -J^, T > T(b). 

In what follows, we shall be working under the assumption obtained by taking the 
conjunction of the assumptions ()4.2()j) and (|4.2fij) . 

After a passage to the standard torus T 2 and an elimination of the xi-variable, we 
obtain an operator 

oo 
3=0 

with 

Po = + e 2 ( S )(0 + ^^ 3 ((t))T(x 2 , + T (e 4 ), 

and 

Pi = e(( qi )) T (x 2 ,0 + O T (e 2 ). 

At this point we may repeat the argument of case 2 to conclude that there exists 
B = & (x 2 , hD x , e, - , T) such that 

^ 0( i )+0 (M +0t ( ( ^ 
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and such that with 



r (x 2 , e, -, T) = (s)(0 + ie((t)) T (x 2j £) + T (e 2 ) + ^(( qi )) T (x 2 , {) + %p 2 {x 2 , £, e, T), 
the symbol 

00 ^-fc 



r oexp (# 6o ) = ^2~^r 



k\ 

k=0 

is independent of x 2 . The construction of the operator B\ is the same as before, 

~ ~2 

and we then introduce a globally defined IR- manifold A C T*T which agrees with 
exp (H bo ) (neigh(£ = 0,T*T 2 )) near £ = 0, is equal to T*T 2 further out, and with 

the property that along A we have Im£i = and 



The action of P e on the corresponding Hilbert space H (A) is microlocally near £ = 
unitarily equivalent to the conjugated operator acting on the L 2 -space of Floquet 
periodic functions. The Weyl symbol of the conjugated operator is independent of 
x and has the form 

p(fi) + e 2 (r (£, e, T) + hn{£, e, T) + . . . 

with 

ro = (s}(0 + ie(((t)) T }(0 + -«<ft»T>(0 + fe 2 + ^ 



e v e 2 

and 

r 3 =0 T {\), j>l. 
The imaginary part of P e along A is given by 

e 3 ((t» T (Rex 2 , Re£) + O (£) + helm (( qi )) T (Rex 2 , Re£) + O {^\ + T (e 4 ). 

Here e ~ /i 1 / 2 . It follows from the assumptions (|4.20J1 and ()4.26|) that choosing T 
large enough but fixed, we can achieve that away from a small neighborhood of 
the torus, the modulus of ImP £ — e 3 ((t)) OQ — /ie((Imq , i)) 0O is bounded from below 
by (e 3 + he)/0(l) ~ e 3 /0(l), when the attention is restricted to the region where 
Re P e is small. At this stage, the situation is completely analogous to the previously 
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analyzed case 2, and the rest of the argument goes through as in that case, without 
any change. 

The following is the main spectral result in case 3. 

Theorem 4.5 Let us keep the general assumptions of Theorem 4.2. In particular, 
assume that the subprincipal symbol ofP e= Q vanishes, so that ifp\ t€ is the subprincipal 
symbol of P e then p l e = eqi + 0(e 2 ), where q% does not depend on e. Let e ~ h 1 ! 2 . 
Assume next that the corrections (s) and (t) are real on the real domain, and that 
the differentials of p and (s) are linearly independent along the set 

Ao,o :p = 0,(s) = 0. 

We assume that Ao,o is connected and that T(0) is the minimal period for the Hp- 
flow in A 0> o- As usual, we shall write S G R 2 and a E Z 2 to denote the classical 
actions and Maslov indices along the fundamental cycles in A , with the first cycle 
corresponding to a closed H p -trajectory. When ((t))oo an d ((Imgx))oo are defined 
as the mean values of (t) and (Im^i) along A 0j o> respectively, we make the global 
assumptions (|4.2U|) and (|4.26|) . Let C > be sufficiently large. Then the eigenvalues 
of P e in the domain 

1 i e 3 
\Rez\ < — , \hnz - e 3 ((t))oo - he((lm <?i))oo| < ^ 

are given, modulo 0(h°°), by the formal quasi- eigenvalues associated with the La- 
grangian torus A ,o, 

p { h (*■ " t) " I) + fl t h ' r > ( h ( k -l)-T^-l)' k= e * ■ 

5 Barrier top resonances: the case of 1 : 1 reso- 
nance 

Consider 

P = -h 2 A + V(x), p(x,0=e + V(x), (x,OeT*R 2 , (5.1) 

where V is an analytic potential, satisfying the same general assumptions as in 
section 7 of ^Hj, which allow us to define the resonances of P in a fixed sector 
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in the fourth quadrant. As in [T^], we assume that V(0) = E > 0, V'(0) = 0, 
V"(0) < 0, and that (0,0) is the only trapped if p -trajectory in p~ 1 (Eq) H R 4 . The 
Taylor expansion of p(x, £ ) in suitable linear symplectic coordinates has the form 

2 A 

p{x, £) - £ = ^ - x ]) + M*) + P4(z) + • • • , - 0. (5.2) 

i=i 

Here Xj > and Pj(x) is a homogeneous polynomial of degree j > 3. We are 
interested in resonances of P near £q, and from we recall that the study of such 
resonances can be reduced to an eigenvalue problem for P — E near 0, after the 
complex scaling given by x = e t7T ^x, £ = e -Mr//4 £, (x,£) G T*R 2 . Performing the 
scaling and dropping the tildes from the notation, we get a new operator with the 
leading symbol 

- M^,0 + te 3 ^ 4 p 3 (x) - tp 4 (x) + ...)=: -g(x,0, (z,£) - 0, 

with 

P2(x,0 = Et^ + ^)' (5 ' 3) 

3=1 

and with the vanishing subprincipal symbol. 

We shall be interested in eigenvalues E of the operator 

Q(x,hD x ;h) =q(x,hD x ) + 0(h 2 ), 

with ~ e 2 , /i* 5 < e <C 1, < 5 < 1/2. After a rescaling a; = ey, we get 

1 1 ~ ~ h 

— Q(x, hD x ; h) = —Q(e(y, hD y ); h), h = — < 1, 

with the corresponding symbol 

-^Q(e(y, r])) ~ ^g(e(y, 77)) + e 2 /?^ (e(y, //)) + .... 

The leading symbol becomes 

-U(e(y, ^)) = P2(y, ^) + iee Zm/i pz{y) - it 2 Pi{y) + ..., 

to be considered in a region where \(y,r])\ ~ 1, where the corresponding eigenfunc- 
tions are concentrated. 

A straightforward application of Theorem 3.1 together with the scaling reduc- 
tions above gives the following result. 
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Theorem 5.1 Assume that the principal symbol p(x , £) of the operator (|5.1|) has an 
expansion (|5.2j) with 

A • k = 0, /or some 0^6Z 2 , (5.4) 

and assume that the average of p 3 along the closed H P2 -trajectories, (ps), vanishes 
identically. Introduce the function 



i r 

s = -pi + — J t{p 3 o exp (tH P2 ),p 3 } dt, 



where T > is the period of the H P2 -flow on p 2 (1), and assume that (s) is not 
identically zero. Then the resonances E of the operator (j5.1|) in the domain 

{z e C;h 2/3 < \z-Eq\ = 0(l)h s }\ [j{z G C; \Rez - E - A \lmz\ 2 \ < 7]\lmz\ 2 }, 

(5-5) 

where 5, r\ > are arbitrary but fixed, are given by 

E = E -ie 2 p(l(k-^-^,e,^ ] h^+O(h™), h=^,keZ 2 . (5.6) 

Here we choose e > with \E — E \ ~ e 2 and the union in ()5.5|) is taken over the 
set of critical values of (s) restricted to p 2 with A varying over this set. 
The function P(£, e,-;h) in ()5.6|) has an asymptotic expansion as h — > 0, 

P(e,6,-;/ i )~p(6) + e 2 ^^k,6, 

i=o V e . 

where _ _ 

r = z(s) + O f e + r, , = O (e+ j > 1. 

In section 7 of [T3] we considered the case of the resonant frequencies (Ai, A 2 ) = 
(1,2), and we have also remarked there that in the case of the 1 : 1 resonance, 
(Ai,A 2 ) = Const(l,l), it is true that (p 3 ) = 0, for any cubic polynomial p 3 . We 
shall now illustrate Theorem 5.1 by discussing this explicit example, where we shall 
take (Ai, A 2 ) = (1, 1). In doing so, we shall also assume for simplicity that p^ = 0. 

It will be convenient to work in the symplectic coordinates (y, rj) given by 



37 



In these coordinates we have p 2 = Y^j=i ^jVjVj: an d the flow is given by 

exp {tH P2 ){ y ,rj) = (e u ^y 1 ,e it ^y 2 ,e^r ]l ,e' UX2 r ]2 ). 
When | a | = 3, we write 

p 3 (x) = X a = ^ a kay k V 



a—k 

J'kaV 'I 

0<k<a 



where 



^kct 



2M/ 2 \ k 
When computing 

1 f 27T 

G = — / tp 3 o exp (tH p ) dt, 
^ Jo 

we use that 

J_ [ 2 \jtx.(2k- a ) dt fTrif A.(2*-a) = 0, 
2vr7 

and obtain 



27r./n ' I i( H 2k- a) ) otherwise, 



0<fc<a 0<fc<a: V J 

2k—a-\-nkQ 2k — a^nkQ 

which we write as 

Go= 9kay k V a ~ k - 

0<k<a 

Here k e Z 2 satisfies (|5.4|) and has the minimal norm and positive first component. 
We are then interested in computing s = (l/2)Ha p 3 , which is equal to 

\h GoP3 = \Y. 9k a a~ ka {y k V a - k , y k , V a ~ k } (5.7) 

0<fc<a 

0<fe<Q 

1 2 ~ 

= 2 Yl 9k a a~ ka a ( k i> <*j ~ k v a i ~ k i)y k+k V° k k ej - 

0<k<a j = l 

0<k<a 

Here a is the symplectic form on R 4 and e\ = (1,0), e 2 = (0, 1). Now when com- 
puting the flow average of (|5.7jl . we notice that the only non- vanishing contribution 
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will come from the terms for which 2(k + k) = 2a + nk for some n € Z. Therefore, 

11 2 

(*> = 3(^3) = 2 E fc^E^'fe - ^)l/ fc+fe - e ^ 2Q - fc - fe -^. (5.8) 

0<fc<a j=l 

0<k<a 
2(k+k) = 2a+nk 

Here we have also used that 

We shall now consider the case when (Ai, A2) = (1, 1) and a = (\a\ , 0) so that 
p 3 (a:) = = 2^ 3/2 (3/? + 3%^! - 3 yi r]l - irfi) . 

Then 

0(3,0)5(3,0) = 2~ 3/2 , a(2,0);(3,0) = 2~ 3/2 3i, 
O(l,0);(3,0) = -2 _3/2 3, O(o,0);(3,0) = -2~ 3/2 i 

and since A • (2k — a) = 2 \k\ — \a\ 7^ 0, we see that 

9ko 



i(2k x - 3) 

Now ko = (1, —1), and it follows that in (|5.8j) we have n = and + A; = a. We get 

(HG0P3) = a ka a a-k,a, 

0<k<a 

and a straightforward computation shows that 

We then see that the differential of this function restricted to the sphere 

vanishes along two closed H p2 -trajectories given by 

xi = Zi = 0, -^7^ = 1, (5.9) 
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and 

*2 = 6 = 0, ^y^ = l, (5.10) 

with the corresponding critical values and 15/4. Notice also that the critical trajec- 
tory (|5.9|) is degenerate in the sense that the transversal Hessian of (s) is degenerate. 
An application of Theorem 5.1 then gives a description of all the resonances of P in 
an energy shell of the form h 2 ^ 3 <C \z — Eq\ = 0(h 5 ), 5 > 0, after we have deleted 
arbitrarily small parabolic neighborhoods of the curves 



{z E C;Rez = E + A |Im z\ 2 }, 



where A G {0, 15/4}. 



6 Complex perturbations on the 2-sphere 

On the sphere S 2 = {x G R ; = 1} C R 3 , equipped with the standard 

metric, let us consider an operator of the form P e = —h 2 A + ieq(x), where q(x) is a 
real- valued analytic function. We introduce the leading symbol of P t =o, p, on T*S 2 , 
and recall that the geodesic flow, which can be identified with the Hamilton flow of 
p, is periodic on each energy surface p~ l (E), E > 0, with the minimal period 

T(E) = -fe (6.1) 

We shall now describe the if p -flow on T*S 2 in a more explicit way. In doing so, let 
us remark that the Euclidean metric on R 3 allows us to embed T*S 2 symplectically 
into T*R 3 , T*S 2 ~ TS 2 ^ TR 3 ~ T*R 3 , so that 

T*S 2 ~ E := {(x, G T*R 3 ; h(x, = 0, h 2 (x, = 0}, 

where hi = x 2 — 1 and h 2 — x ■ ^, and where we use the Euclidean scalar product 
and the norm. In what follows we shall write {•, ■} to denote the Poisson bracket on 
T*R 3 and the Poisson bracket on S will be denoted by {■, 

S is a symplectic submanifold of T*R 3 , since {hi, h 2 } = —2 there. The restriction 
to S of the Liouville form £ ■ dx can be identified with the corresponding Liouville 
form on T*S 2 . We then have the corresponding identification of the symplectic 
forms. Furthermore, p = £ 2 , when viewed as a function on S, and the corresponding 
Hamilton field satisfies Hp = H v mod TS CT , where H p is the Hamilton field on T*R 3 
and the exponent a indicates that we take the symplectic orthogonal. It follows that 
Hp = Hp, where p G C°°(R 3 ) is such that p = p along £ and Hp is tangent to E. 
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Writing p = p + ah\ + bh 2 , where a and b are chosen so that {p, hj} = along E, 
j = 1, 2, we find after a simple computation that p = x 2 ^ 2 . Restricting to x 2 = 1, 
we get 

-H x 2 e = £ ■ d x - i 2 x ■ % 
This implies that the map exp {tH^/2) : T*S 2 \0 -> T*S 2 \0 is given by 

(x, ^ (cos(|e| t)z + sin(|£| f) j|, - |£| sin(|£| t)z + cos(|£| *)£)■ (6.2) 

Example. When q(x) = X1X2, we shall compute the averaged perturbation 

{q)=i LL v° ex v( tH ?/ 2 ) dt ° n p' 1 ^' ( 6 - 3 ) 



and using ([6.2)1 we obtain 



(?>(af,0 = o • ( 6 - 4 ) 



Now (g) is invariant under the if^-flow, and should therefore be viewed as a function 
on the reduced space of oriented closed orbits in 

0:=p-\l)/exp(RHZ). (6.5) 

Here O is a two-dimensional compact symplectic manifold, which consists of all 
oriented great circles of S 2 , and we can identify O with S 2 by associating the unit 
vector y = x x £ to the Hlf— trajectory in p _1 (l), given by ([6.2)1 . In order to express 
the flow average (q) in terms of y, we notice that we can choose (x,£) corresponding 
to y with 

= H/2,yi,o) (6 6) 

(2/f + 2/l) 1 / 2 ' 1 J 

and 

It follows that 

(g>(y) = — ■ (6-8) 

A straightforward computation shows that (g) is a Morse function on S* 2 , with 6 
critical points. We get two non-degenerate maxima at ±(l/ v / 2,-l/v / 2,0), with the 



41 



corresponding critical value 1/4, two non-degenerate minima at ±(1/V2,1/V2,0) 
with the critical value —1/4, and two non-degenerate saddle points at the poles 
±(0, 0, 1), with the corresponding critical value 0. 

Remark. The considered example comes essentially from Notice also that the 
expression (|fi.8j) is in agreement with the general fact, established in jHj, that the 
Radon transform q \— > (q) on S 2 maps the space of the restrictions of the homoge- 
neous harmonic polynomials of a fixed degree into itself, and is a multiple of the 
identity when restricted to this space. A simple proof of this fact following [Hj, is 
obtained if we observe that the Radon transform commutes with the SO (3) action 
on S 2 , and apply Schur's lemma. Let us also recall from ^T] and [H] that the kernel 
of the Radon transformation, viewed as a mapping on C°°(S 2 ), consists precisely of 
all odd functions on S 2 , and that it becomes bijective when restricted to the space 
of smooth even functions. It follows that every smooth Morse function in the range 
of the Radon transform has at least two saddle points. 

As a preparation for further considerations, we shall now derive an expression for 
the Poisson bracket of two functions f,g G C°°(T*R 3 ), viewed as functions on the 
symplectic manifold £. In doing so, we have to compute the Hamilton field of the 
restriction of / to £, and repeating the previous arguments, we find that Hj = Hj, 
where 

{h 2 ,hi} \hi,h 2 \ 



Therefore, 



{/, 9h = {/, 9} + \ ({/, h 2 }{hi, g} - {/, h}{h 2 , g}) . (6.9) 



We summarize the discussion above in the following essentially well-known pro- 
position — see also jl] for a more general discussion of Hamilton mechanics with 
constraints. 

Proposition 6.1 Let 

£ = G T*R 3 ;fci(x,0 = 0,/i 2 (x,«£) = 0}, hi = x 2 - 1, h 2 = x ■ f. 

Then £ is a symplectic submanifold o/T*R 3 and the Hamilton vector field ofp = £ 2 , 
viewed as a function on £, is given by H x 2^2. The Poisson bracket of the restrictions 
off,ge C°°(T*R 3 ) to £ is given by (loTIJ) . 

We now consider the case when the perturbation q is an odd function in R 3 . Then 
(q) = 0, and assuming that q is an odd monomial, we shall derive an expression for 
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the second averaged correction (s), given by (|2.6jl . In doing so, when (x,£) G T*S 2 \0, 

l p/2 



we introduce z = x + i£/ |£| G C 3 . It follows then from (|6.2|) that along the H^ 



trajectory we have 

z (t) = e-'^(0). (6.10) 

Therefore, with q(x) = x a , \a\ odd, we get 

«K«q,(t/2flf = ( ze "" s " 2 + ^'T 4S(;) ,wi«h-iw. 

We now have to compute the function Go, given by ()2.2j) . 

G = j^j tq o exp (tH p ) dt on 
and we immediately see that 

G ° = 9N+1 m J2 ( < n) ^uiloiAiv (6 - n) 



We shall next consider the Poisson bracket of Gn and 



9 M / ' 



a 



so that 



0<7<a V 



0</3<a 
0<7<Q 



When computing the flow average of this expression, we see using (jfi.lOJ) and the 
fact that |£| is constant along the Hp -Raw, that the only non-zero contributions to 
the average come from the terms for which |a| = \/3\ + |7|. We get 



0</3<q 
0<7<a 
*I = I/3| + M 



We shall now illustrate the preceding discussion by an explicit example when 
a = (1, 0, 0), q(x) = X\. Using (|6.11|) we see that 



2£ 



2 ' 
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and we may then check directly that H^Gq = q. A direct computation using ()6.9|) 
shows next that 

i^o,g|s- - y—^ + -^r 

When computing the flow average of {Go,q}s restricted to we notice that 

x\ + £f is flow invariant, while using (J6.2j) we find that the flow average of £ 2 is 
{x\ + CD /2. Therefore, 

<s> = -^<{Go,g}£> = ^-^^, on p-i(l), 

and using ()6.6|) and (j6.7|) . we find that when viewed as a function on O ~ S* 2 , (s) 
becomes 

<*> = i/i 2 -^ (6-12) 

The differential of this function vanishes along the equator y\ + y\ = 1, with the 
corresponding critical value —1/8, and at the poles ±(1, 0, 0) with the critical value 
1/4. ' 

Combining Proposition 2.1 and remark at the end of section 2, together with 
Theorem 4.2 and the discussion above, we obtain the following result. 

Theorem 6.2 When q(x) = X\ on R 3 , let us consider the operator P e = —h 2 A + 
ieq(x) acting on L 2 (S 2 ). Assume that e <C h 1 ^ 2 . Then the spectrum of P e in 
neigh(l, C) is contained in the union of the rectangles of the form 

h 2 (k + + [~0(e 2 + h 2 ), 0(e 2 + h 2 )} + i[-0(eh), 0(eh)}, k e Z. 

Let F G (—1/8,1/4) and restrict the range of e further by assuming that /i < e. 
Then the eigenvalues of P e in the set 



Rez - h 2 ( k + ^] - e 2 F 



are given by 



j=0 
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Here 

roU,e,~^ = (sKZ) + o(e+±Y r, U,e,^j = O (e+^j , j > 1, 

are holomorphic in £ G neigh(0, C 2 ) and smooth in e,h/e G neigh(0,R). The coor- 
dinates £1 = ^i(-E) and £ 2 = ^{E, F) are the normalized actions of 

A E , F :p = E,(s) = F, 

for E G neigh(l, TV), F £ neigh(F , R), given by 

= 77- I / vdy- vdy), j = 1,2, 

with jj(E,F) being fundamental cycles in Ae,f, such that ji(E,F) is given by a 
closed Hp-trajectory of minimal period T(E). We have £i(E) = E 1 / 2 — 1, so that 
when expressed in terms o/£i, p becomes (£1 + l) 2 . Finally, a G Z and S G R are 
fixed. 

Remark. We shall finish this section by pointing out that the methods of and of 
the present paper can be extended to cover the case when the operator P e , introduced 
in section 1, depends holomorphically on a parameter z G neigh(0, C), in such a way 
that on a symbol level, we have a holomorphic family of functions P(x, £, e, z\ h), with 
all the properties described in section 1 holding uniformly in z, and such that the 
leading symbol of P(x, £, 0, z\ h) is of the form p(x, £) — z, for a real- valued p(x, £) 
with a periodic Hamilton flow. We then check by inspection that the reductions 
of sections 3 and 4 work in the same way, and as a microlocal Birkhoff normal 
form for P e , we obtain a family of translation invariant operators on T 2 depending 
smoothly on z G neigh(0, C). The setup of the global Grushin problem is the same 
as before — see also |IBj, where the general holomorphic dependence on the spectral 
parameter is assumed from the beginning. The complete asymptotic expansions 
for the eigenvalues then come from applying the implicit function theorem to the 
"effective Hamiltonian" E_ + , as in section 6 of [THj . 

These observations are motivated by the problem of studying asymptotics of 
eigenfrequencies associated to the damped wave equation with an analytic damping 
coefficient, on an analytic Zoll surface M. Indeed, let us recall from [201 anc ^ 
that after a semiclassical reduction in the eigenfrequency equation 

(-A + 2ia(x)r - r 2 ) u = 0, |r| > 1, 



45 



obtained by writing 



r=-V, 0</i<l, 2! G neighfl, C), 

n, 

one is led to consider the problem 

(-h 2 A + 2ih^fza(x) - z) u = 0. 

Therefore, in this case P e = —h 2 A + Ht^fza with e = h, and the subprincipal 
symbol of P e= o = —h 2 A vanishes. Without any analyticity assumptions and without 
any restriction on the dimension, it was proved in [Hj that the spectrum of P e in 
neigh(l, C) has a cluster structure, for any a G C°°(M; R). If now M is an analytic 
compact symmetric surface of rank one which is not S 2 , then from JU] we know 
that the flow average (a) ^ unless a vanishes identically. The results of [To] , 
generalized to cover the case of holomorphic dependence on the spectral parameter in 
the non-selfadjoint perturbation, as indicated above, give then complete asymptotic 
expansions for the eigenvalues z, in the sub-clusters corresponding to regular values 
of (a), viewed as a function on and to non-degenerate extreme values of 

(a) on p -1 (l). If M = S 2 and a is not an odd function, then the results of ^3] 
are still applicable, since, as was remarked before, in this case (a) does not vanish 
identically. Consider finally the case when M = S 2 and a is an analytic odd function. 
It follows then that the results of sections 4 and 6 of the present paper, generalized 
as described above, give complete asymptotic expansions for the eigenvalues z in the 
sub-clusters corresponding to regular values of the averaged second correction (s), 
after the damping coefficient a has been multiplied by a sufficiently large but fixed 
coupling constant. 
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